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Abstract. We consider a nonlinear evolution problem with an asymptotic parameter and construct 
examples in which the linearized operator has spectrum uniformly bounded away from Re z > 
(that is, the problem is spectrally stable), yet the nonlinear evolution blows up in short times for 
arbitrarily small initial data. 

We interpret the results in terms of semiclassical pseudospectrum of the linearized operator: 
despite having the spectrum in Re z < -yo < 0, the resolvent of the linearized operator grows very 
quickly in parts of the region Re z > 0. We also illustrate the results numerically. 



For a large class of nonlinear evolutions the size of the resolvent has been proposed as an 
explanation of instability for spectrally stable problems. Celebrated examples include the plane 
Couette flow, plane Poiseuille flow and plane flow - see Trefethen-Embree |l6l Chapter 20] for 
discussion and references. Motivated by this we consider the mathematical question of evolution 
involving a small parameter h (in fluid dynamics problem we can think of h as the reciprocal of 
the Reynolds number) in which the linearized operator has the spectrum lying in Re z < -yo < 0, 
uniformly in h, yet the the solutions of the nonlinear equation blow up at time 0(1) for data of 
size 0(exp(-c//z). 

We know of one rigorous example of such a phenomenon given by Sandstede-Scheel fT2]. 
They considered Ut — Uxx + Ux + u^ on [0, i] with Dirichlet boundary conditions, and showed that 
blow up occurs with arbitrarily small initial data as ^ — > oo. In that problem h - The paper 
f|l2i is our starting point and we use its maximum principle approach to obtain results for suitable 
operators in any dimension. In addition, we emphasize the connection with the semiclassical 
pseudospectrum and provide some numerical comparisons. 

We consider a semiclassical nonlinear evolution equation 



1. Introduction 



(1.1) hilt = P{x,hD)ii + 11^, xeR'^ , t>0. 

where P{x, hD) is the following semiclassical differential operator 
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for some C, Ca,k,M > 0. Also, p £ C°° has the properties 



(1.4) l^pl |Vp|2>4(^ + ,.„), 



for some Cij,yo,N > 0. Finally, ^> 0. 

Remark: All of our results hold for weaker assumptions on the growth of V and p, however ( |1.3| l 
and ( |1.4| l are convenient for our purposes. 

We will show in section |2] that for V{x) and p as in (1.3 1 and ( |1.4| respectively, the linearized 



problem is spectrally stable, that is, the spectrum is bounded away from Re z > uniformly in h. 
Yet, we also show that h as an unstable equilibrium at u = for all potentials V{x) satisfying 
(|1.3| and all p satisfying (|1.4|. Specifically, we show 



Theorem 1. Fzz /i > 0. Then, for each 

<h <ho, 

where ho is small enough, there exists 

meC^iU"), uo>0, INIb <exp(-^), p = 0,l,..., 

such that the solution to ^1.1) with u{x,0) = uo{x), satisfies 

— > oo, t — >T, 

where 

T = 0(1). 

A nice example for which our assumptions hold is \l.2\ with x e R, V{x) — x^, and ( Vp, D) - Dx- 
That is 

(1.5) Piix, hD):=- ({hDxf + x^) + ihD^ + ji, xeR 

It is easy to see (and will be described in Section |2]) that 

Spec(Pi(x,/jD)) = - 1/4 - hiln + 1) : n = 0, 1,2 • • • } 
c {z : Rez < /i - 1/4}. 

For n > \ the spectrum intersects the right half plane and thus instability of the linear problem 
follows. We are interested in the range < /i < |, where we will relate the instability of u = to 
the presence of pseudospectrum in the right half plane. 

For more about ( |1.5| l see p6'. Chapter 12]. In particular, Cossu-Chomaz [Ij relate it to the lin- 
earized Ginzburg-Landau equation and analyze resolvent of ( |1.5| l and the norm of the semigroup 
gPi{x,hD)t numerically. 

The operator ( |1.5| l is also closely related to the advection-diffusion operators mentioned above, 
—Dy + iDy = dy + dy, on [0,/"], with, say, Dirichlet boundary conditions; see fS", Chapter 12] for 
a discussion and references. When rescaled using x - yjl, h = l/l the operator becomes the 
semiclassical operator -QiD)^ + ihDx on [0, 1]. When the domain is extended to R, the potential 
x^ is added to (hDx)^ to produce a confinement similar to a boundary. 

We relate the blow-up of solutions to ( |1.1| | to the presence of pseudospectrum of ( |1.2| | in the 
right half plane. However, because estimates on semigroups for ( |1.1[ | with quasimode initial data 
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are poor, we are unable to exhibit blow-up starting from a quasimode. Instead, we present a 
simple and explicit construction of quasimodes for P{x, hD) (for a more general setting see IISll. 
We then use these quasimodes as initial data in numerical simulations and observe that, although 
in some cases the ansatz solution blows up more quickly, the solutions with quasimode initial 
data behave similarly to what is expected from a pure eigenvalue for ( |1.2| | with positive real part. 




Figure 1. The plot shows numerical simulations of the evolution of with 
h - 1/193 and two initial data. The evolution with initial data a real valued 0{h'^) 
error quasimode with eigenvalue z = is shown in the top two graphs and that 
with the ansatz constructed in the proof of Lemma |3] as initial data is shown in 
the bottom two graphs. We observe that, when the initial data is a quasimode, 
blowup occurs in time « 0.3, while for ansatz initial data, blow-up occurs in time 
« 0.175. However, as would be expected from eigenfunction initial data, we see 
that the solution with quasimode initial data exhibits little transport to the left. 
On the other hand, the ansatz transports left significantly. 

The paper is organized as follows. In Section |2] we review the definitions of spectra and 
pseudospectra and discuss them for our class of operators. In Section |3] we give a construction 
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of quasimodes for one dimensional problems. Although the results are known, (see ||3|,||5|, ||T31 ) 
a self-contained presentation is useful since we need the quasimodes for our numerical exper- 
iments. Also, there is no reference in which analytic potentials (for which quasimodes have 
0{exp{-c/h)) accuracy) is treated by elementary methods in one dimension. Section|4]is devoted 
to the proof of Theorem [T] using heat equation methods. Finally, in Section|5]we report on some 
numerical experiments which suggest that quasimode initial data gives more natural blow-up 
and that blow-up occurs at complex energies. 

AcKNOwLEDGEMNTs. The author would like to thank Maciej Zworski for suggesting the problem 
and for valuable discussion, guidance, and advice. Thanks also to Laurent Demanet and Trever 
Potter for allowimg him to use their MATLAB codes, Justin Holmer for informing him of the 
paper by Sandstede and Scheel, and Hung Tran for comments on the maximum principle in 
Lemma |3j The author is grateful to the National Science Foundation for partial support under 
grant DMS-0654436 and under the National Science Foundation Graduate Research Fellowship 
Grant No. DGE 1106400. 



2. Spectrum and Pseudospectrum 



We do not use the results of this section to prove Theorem T] Instead, we present them to 
emphasize the connection of the size of the resolvent with instability. We believe that instability 
based on quasimodes would be more natural and allow for proof of instability at complex 
energies. We illustrate this with numerics in Section |5j 

To describe the spectrum of P{x, hD), we observe that 



e 2h P(x, hD)e 21, - 



(/2D)2 + y(x) + -|Vp|' 



Thus, the spectrum of P{x, hD) is given by that of a Schrodi nger operator with potential V{x) + 
||Vpp -I- |Ap. Since V{x) and p have the properties given in |l.3l and (1.4 1 respectively, P{x,hD) 



has a discrete spectrum that has real part bounded above by -70 (see for instance |8l Section 6.3]) 



We now examine the pseudospectral properties of ( 1.2 1. 
Definition. Let Q{x, hD) be a second order semiclassical differential operator. Then, z e A(Q) if 
and only if 3 u{h) £ H^iW^) such that ||w||l2 = 1 and 

\\{Q{x,hD)-z)u{h)\y=0{hn. 



We say z is in the semiclassical pseudospectrum of Q if z e A(Q). 

Remark. We note that for z € A(Q), \\Qix,hD) - z)-i|| > h-^/CN, for any N. This relates 
our definition to the more standard defintions of pseudospectra in terms of the resolvent. For 
discussion and generalizations see Dencker [4J and Pravda-Starov lUTl . 

The criterion for z £ A(Q) is based on Hormander's bracket condition (see Zworski IITSl and 
Dencker-Sjostrand-Zworski Il5l): 

(2.1) Q{xo, ^0) = z and {Re Q, Im Q}(xo, 4o) < 0, 

then z e A(Q). We use this condition to show that the pseudospectrum of P{x,hD) nontrivially 
interesects the right half plane. Specifically, 
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Lemma 1. For P{x, hD) given by (|L2|, A{P{x, hD)) Clilmz = 0} = (-co, 

Proof. First, observe that 

P(x, E) = -lEl" + zXVp, O - Vix) + ^ and {Re P, Im P} = -2(5^-^' + (VV, Vp). 

We have assumed Imz = 0. Therefore, we need only show that, for a dense subset U c (-co, ^], 
y e U implies that there exists x such that ( |2.1| holds for the symbol P(x, (^), at (x, 0) with z - y. 

We proceed by contradiction. Suppose there is no such U. Then, there exists O c [0, oo) open 
such that for all x e V'^iO), (VV, Vp)(x,0) > 0. Let cpt := exp(fz<Vp,D)) be the integral flow of 
f<Vp,D) and £ R'*havey(xo) = 0. Define /(f) := V{(pt{xo)). Then(9f/ = <Vy((pt(xo)), Vp((pt(xo))). 



Suppose that (pt{xo) escapes every compact set as \ t\ increases. Then ( |1.4[ ) implies that f{t) — > oo 
as \t\ increases. Let iv e O and to ini{t e R : f{t) - w]. Then to is finite since w > /(O) and 
f(t) — > oo. Together, f(to) = w e O and /~^(0) open imply the existance of 6 > such that for 
te{to- 6, to + b), fit) € O. But, f{t) e O implies f'{t) > 0. Therefore, f{t) < w for t e {to - b, to) 
and thus, since f{t) — > oo, there exists t <to such that f{t) -w,a contradiction. 

We have shown that there is a dense subset U c (-oo, y\ with U c A(P). Hence (-oo, c A(P). 
Next, observe that sup Re P(x, (^) = ^lj and thus, A(P(x, /zD)) Plllin z = 0} = (-oo, ^\ as desired. 

To finish the proof, we need only show that (pt{xo) escapes every compact set. Suppose the 
flow at xo exists for all f 6 R. Define h{t) := p{(pt{xo)). Then dth = |Vpp > c > and we have that 
h — > +00 as f — > +00. But, p £ C°° and is therefore bounded on every compact set. Thus, (pt{xo) 
escapes every compact set as f ^ +oo. 

Now, suppose the flow at xo is not global. Then, <pf (xq) is an integral curve of /(Vp, D) with t 
domain a proper subset of IR. Thus, as proved in [lOl Lemma 17.10], (pt{xo) escapes every compact 
set. □ 

Putting this together with our discussion of the spectrum of P{x, hD), we have that for < /.z 



and p as in \\A) , although Spec(P) is bounded away from Re z > 0, A(P) nontrivially intersects 
Rez > 0. 

For the specific case, V{x) - and Vp constant with |Vp| = 1, the above argument gives us 
that 

spec (P(x, hD)) = |-(2n + l)/2 + (^ - ^) : n > o| . 
In addition, the pseudospectrum is given by, 

A{P{x,hD)) = jz : Rez < -(Imz)^ + ^) . 

We see that for ii> \, the spectrum interesects the right half plane and so instability of w = is 
a classical result. However, for < /j < |, the spectrum is bounded away from the Rez > and 
only the pseudospectrum enters the right half plane. Yet, in the regime < /j < |, we will show 
that M = is unstable and, moreover, for arbitrarily small initial data, the solution blows up in 
finite time. 
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H-1/4 3z 




Figure 2. We see that the spectrum of \1.5\ (blue dots) is bounded away from 
Re z = 0, while the pseudospectrum (shaded region) enters the right half plane. 
The region for which we prove blow-up corresponds to the dashed red line. 



3. One Dimensional Quasimodes 

We proceed by constructing quasimodes for operators in the one dimensional case with 
z( Vp, D) = dx- We implement WKB expansion for the quasimode following the method used 
in in. Let 

(3.1) P{x,hD) - (hDxf + ihD^ + V, 

where V € C°° and V may be complex. 

Remark. The following theorem is a special case of general theorems about quasimodes flSj 
Theorems 2 and 2']. For the reader's convenience we present a direct proof in the spirit of Da vies 
0. 

Theorem 2. Suppose that P{x, hD) is given by ( |3.1| l and that 

rl 



z = -^l + i^o + V{xQ), 



where xq satisfies the condition that 



Rey'(xo) > 0. 

There exists an h-dependent function cp e C^(R), such that \\q)\\i2 = 1 and 

mx,hD)-z)cp\y = Oih-'). 

In addition cp is microlocalized to {xq, <^o) in the sense that for every g e C'^{W x ]R") vanishing in a 
neighbourhood of{xo, <^o)/ 

\\g{x,hD)(p\y = o{hn. 
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When V is real analytic than we can find (p such that 

\mx,hD) - z)(plb < Cexp(-l/C/i) . 

Proof. Let x ^ C^(]R) with xi^) = 1 if 1^1 < 6/2 and x{x) = if |x| > 6 where 6 will be determined 
below. Define / := exp{ixp/h)a{x) where 

N-2 

a{x) = ^ am{x)h^. 

m=0 

Finally, let g{xo + x) := x{x)f{x) for all x e R. 

We will find appropriate am and xjj in what follows. First, by a simple computation 

/■ N \ 



{P{x,hD)-z)f = 



"'cb 



where <pm are inductively defined by 

4)m := {-{il^'f + if + V- z)am + if'am-i + (2zi/;' + '^)a'ni-i + C-v 

where we use the convention that fl^ = for m > N - 2 or m < 0. Now, we set (^^ = for 
< m < N - 1. Given that 6 is small enough, this will enable us to determine all am as well as i/^. 

Observe that, using the condition, (pQ = 0, we obtain 

l/;'2 _ i^jj' = V-Z. 

Now, letting z = -S,^ + i£,o + V{xo), we have a complex eikonal equation 

• 2 • 2 

(^-'2) =^(^o + x)-y(xo) + (4o-0 . 

Then, letting ip - ip - ^x, we have 

xf>{x) := £^V{xo + t)-V{xo) + {^o-'2)^ dt 

= (^0 - + ^ (^0 - ' V'{xo)x^ + 0(x3)j 

and hence 

4(^0 + 1) 

Now, we have assumed that Re V'{xo) > 0. Therefore there exists 7 > such that 

yx^ < lmxl){x) < 3yx^ 
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for all small enough x and h. Also, for x and h small enough 

satisifies \6{x)\ < jS. We choose 6 > small enough so that these conditions both hold for <h< 6^ 
and \x\ < 6. 

The condition <pm+i = 0, implies 

«m = -Q{i^"cim + 

with the convention that = and initial conditions, 

flo(O) = 1, flm(O) = 0, m > 0. 
Putting G{x) := i^p" {y)d{y)dy we obtain aq = exp(-G(x)) and 

(3.2) := -e-GW f ^ > q. 

Jo 

Before proceeding to show exponential error for V analytic, we show 0{h!^) error for arbitrary 
V. To complete the proof of OQ^) quasimodes, we need to estimate 

\mx,hD)-z)g\\i2i\\g\h. 

Let C denote various positive constants that are independent of h and x. Then, 

X6/2 r^l2 
\f{xfdx> e-'y^'^'-'-'^dx 
6/2 J-5/2 

XSr 1/2/2 r-1/2 
^-6yf-C^ll2^i > e-'y''-^h'/^dt = Ch^l\ 

6/1-1/2/2 J-1/2 

Next, we compute 

\mx,hD)-z)gh2 = Wh'fx" +2hYx' + hfx' +x{Pix,hD)f-zf)h2 

(3.4) < h^\\fx"\y + 2h^\\f'x'\y + /j|I//IIl2 + \\h''x<pNe''f"\2- 

Thus, we need to estimate each of the norms. Note that x' and x" have support in {x : 6/2 < \x\ < 
6}. Thus, we have 

(3.5) \\fx"\\l < Cs r e-^y'^^~'+^dx<Ce-y^^'^. 

^ J5/2<\x\<6 

Similarly, 

(3.6) \\f'x% < Ce-y''!"^, \\fx%^ < Ce-y'""^. 
Next, observe that 

g-2y.2ri+c^^ < C/z2N||(^nI|L e-^y-'^'^h^Hx 
e-^y'^^^h^lHx < c(//~||(/)nIIl»P'/'. 

oo 
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Now, \(pm\ < Cm on {x : \x\ < 6}, uniformly for h < 6^. Therefore, combining ( |3.4| | with inequalities 
( |3.3| l, (3.5l, (3.6l, and ( |3.7| |, gives OQi^) quasimodes for arbitrary N. We then normalize to obtain 
(p. 

We will now assume that V{x) is real analytic and prove exponential smallness of the error. 

Lemma 2. Let I = [-6,6] where 6 is a small constant. Suppose that tq, ti, T2, flnd do fl?'^ holomorphic 
functions ofz e Q and |t2| > ^/o'' some C > wfere Q fs a neighborhood of I in C. 
if d„i is defined inductively by 

(3.8) d„Mz)= f T2{Qd':,{Q + TiiQd'JQ + ToiQdUQdy. 

Jo 

Then for some Ci > 0, C2 > and [-6,6] c Q c Q, 

(3.9) sup \dPdm\ < V'C'r^m'". 

Proof. Using integration by parts, we obtain that 

£f„+l(z) = T2{z)d'Jz) - T2(0)C(0) + (T1(Z) - T^(z))d„(z) + f (t^'(C) - T;(C) + To(C))d,„(C)rfC. 

Jo 

Then, since T2 is holomorphic in Q and info |t2| > ^ for some C > 0, we make the conformal 
change of variables z — > w where 



dw 1 



Then, letting bm = d„-i{ziio)), we have 



o'f„(&m+l)(w) = dC bm{zv) - 6po(o'r„&m)(0) + o'f„(po&m)(w) + olf„ I pi{Qbm{Qd(, 



where po = (ti - t'^{z{w)) and pi = ((T2 - + i^o)t2)(z(^))- Put Qa, := {w : z{w) e Q}. Then, since 
the change of variables was conformal, and Ti, i = 0, 1,2, are holomorphic, we have that there 
exists Cp > such that 

|^f„p,b„<CPV for z = 0,1 
where we define |/|q := sup^ |/| for a function / defined on Q. 
We claim that for some Ci > 0, Co > Cp, 

\dlbJn,.,<Cl^'c^^Hm+pr^P. 

We prove the claim by induction. The holomorphy of bo gives us the base case. We now prove 
the inductive case. 

By the inductive hypothesis, we have that 

(3.10) IdC'KK. < Cl^^C^*\m + p + l)'«+P+i. 
Similarly, 

(3.11) |fem(0)h,„ < Cf+^m-. 
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Next, we prove similar estimates for d^{pobrn)- By Leibniz rule, we have that 



(3.12) 
where 



\dnpobm)\o^ = 



r 

E 

k=0 



k\ip-ky. 



— / , 1 Knt,p 



-k\m + p-kr+P-K 



We claim that for < fc < |, rk^m,p ^ Tp-k,m,p- To see this, we write this inequality as 

]^im+p- kr+P-^ >{p- kf-\m + kf+^, for < A: < | . 

Putting x:= ^ and y = the inequality is equivalent to 

x\l + yf+y > 1/5/(1 + x)i+^ for < X < 1/ 

which follows from the monotonicity of the function x (^) (1 + 
Next, observe that, <k < p - 1, 

rk+i,m,p p-k{k + 1)*^+! {m+p-k- i)"'+P-'^-i 



^k,m,p 



k+1 k^ {m+p- kf+P-k 



p - k 



1 V kl \ m + z 



m + p -k-l\ kl \ m + p-k 

p-k 1 1 



m+p-k 



m+p -k-1 

where we use log(l - x) <-x + ^. Then, since for < A; < |, r!^^m,p ^ ^p-k,m,pi we have 

?+i k-i 



fc=0 
H+l 



fc=0 



«=0 



E+1 



k=0 k=0 

< d'^^C^+^ip + 2)rQ,m,pe^ < e^Cl^^C^'-^{m + p + 1)'"+?+^ 



Therefore, there exists Mi > such that 

(3.13) \dl{pMci^ < MiCl^^C^^\m + p + Xf^^^K 

By analagous argument, there exists M2 > such that 



(3.14) 



n-z(w) 

Jo 



.(C)rfC 



dl-\ipid^z)b„,) < M2Cl^^C^^\m + pr^P. 
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Next, choose Ci > Q (4max(Mi,M2, 1)) . Then, combining ( |3l0| , (|3lT), ( [3l3| , and ( |3l4| , we 
have 

|5f„&.+l(lf^)lQ. < C^+Vr+2(^ + p + i)m+p+l_ 

Then, since zf — > z is a change of variables independent of m which maps Cl-a, — > Q and h„i{w) - 
d„,{z{w)), we have 

Now, choose Q c Q with inf{|z - Q : z e Cl, C, e dQ.} > y > 0. Then, since cf„j are holomorphic, we 
apply Cauchy estimates to obtain 

m+l...m 



□ 



We apply the lemma with d„,(x) := £'G(^)fl^(x), T2 := -0, ti := 20G', and tq := 0(G" - (G')^) 
where G and G are given above. Analyticity of V implies that uq, 6, and are holomorphic in a 
neighbourhood of I. 

Then, putting 1/N — eCih, using Lemma|2]and that i/^ is real analytic, we have 



(3.15) 



sup < h^C{C2C^iN-lf-^+2ClC^iN-l)^-^ + 6ClC^{N-2f-^) 

xe[-6,6] 

/ 1 \^ 1 



where C denotes various postive constants that are independent of N. Finally, combining l\3A) 

with inequalities <\3.3) , <\3.5) , <\3.6) , <\3.7) , and ( 3.15| l gives 0(e~cs) quasimodes. We then normalize 
g to obtain cp. □ 



Now, applying this result to Pi{x,hD), we obtain 



2xot + 



\l/2 



for xq <0 and z = -(^^ + iE.o + f-i- x^. 



4. Instability 

Our approach to obtaining blow-up of will follow that used by Sandstede and Scheel in 
||12^. We will first demonstrate that, from small initial data, we obtain a solution that is > 1 on a 
deformed ball in time h = 0(1). We will then use the fact that the solution is > 1 on this region 
to demonstrate that after an additional t2 = 0{h), the solution to the equation blows up. 

First, we prove that there exists initial data so that the solution to is > 1 in time 0(1). 
Recall that q)t denotes the flow of /( Vp, D). 
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Figure 3. We set h = lO"'^ and see that the difference between the solution to 
with initial data a quasimode with error 0{h^) (red line) and the solution with 
initial data a quasimode with error 0(/2'^) (blue dots) is negligible. Thus, by using 
0(h^) error quasimodes, we have not introduced large error into our numerical 
calculations. 



Lemma 3. Fix ^ > 0, a < < e < 5 (/J - a), and {xQ,a,b) e R'* x R x such that both 
(pf{B (xo,2a)) c V~^[0,^ - a - e]for t < 6 and cpt is defined on B {xo,2a)for <t <26. Then, for each 

<h <ho 

where ho is small enough, there exists 

Uo{x)>0, llwollc^' ^ exp^-^j, p = 0,l,... 



and < ti < 6 so that the solution to {l.D with initial data uq satisfies u{x,ti) >1 on x e cpti{B{xo,a)). 



Proof. Let v solve 
(4.1) 



(hdt - Pix,hD))v = 0,vix,0) = VQ. 



NONLINEAR INSTABILITY IN A SEMICLASSICAL PROBLEM 



13 




Figure 4. We show a numerical simulation u{t) of the evolution of with 
a quasimode at imaginary energy as initial data. Specifically, we set h = 1/193 
and use a quasimode with eigenvalue z = + The real part is shown in 
the top graph and the imaginary part in the bottom graph. We see that although 
subsolution methods do not apply to these quasimondes, blow up still occurs in 
time « 0.35. 

Let w : ~^ 3^ and define O :- {x : wq > 0}. We make the following assumptions on wq, 
(4.2) wo>0, llu^olb < exp(--^), wq e C{R'^) 



(4.3) Wo e C°°(0), supp zvq c B{xQ,2a), Wq > exp (^-^j on B{xo,a), 

(4.4) do is smooth, -Awo{x) < Cwo{x) - fiforx e O and <h <ho. 

where C°°(0) are smoothly extendible functions on O. We will construct such a function at the 
end of the proof. 
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Define w : [0, 26) x R'* ^ R by 

iexp (^ft^ zvQ{(pt{x)) where cpt is defined 
else. 

Since supp w c B{0,2a) and (pt is defined on B{0,2a) x [0,26), w is continuous. We proceed by 
showing that is a viscosity subsolution of ( |4.1| in the sense of Crandall, Ishi, and Lions ||2| . 

First, we show that w is a subsolution on Of := (ptiO) for t < 6. 

hivt - P{x, hD)w = hwf - h^Aw - ih{yp, D)w - lav + V{x)w 
= (a - n + V{x))w - h^Aiv 

< exp ^^f j ((a - + V{x))zvo) - h^Azv 

Now, by Taylor's formula, (pt{x) = x + 0{t). Hence -A[wo{q)t{x))] = -Awo{x) + 0(f). We have 
t < 6, and -Azvq < Cwq - |S on O. Therefore, for 6 small enough, -Aw < Cwq. Hence, for h small 
enough independent of < 6 < 6o/ 

hzvt - P{x, hD)zv < expl^^tj (a- [i + Op- + V{x)^ Wq 

Now, since for some e > and t < b, supp w <z V~\Q, - a - e\ we have that 

hwt - P{x, hD)w < exp {^t^ {Ch^ -e)zv <0 

for h small enough. Thus, if is a subsolution on Ot for t < 6. Next, we observe that on (R** \ Of), 
w = and hence is a subsolution of ( |4.1| |. 

Finally, we consider ^Of := (pt{dO). We have that ^Of is smooth. If i/o ^ dOt and w is twice 
differentiable at yo, then Wt - {Aw){yo) - {Dw){yo) - w{yo) - and w is clearly a subsolution to 
( |4.1[ | at yo- Let yo € dOf be a point where w is not twice differentiable. Suppose that (p e such 
that w — (f) has a maximum at i/o- 

We take paths through yo to reduce to a one dimensional problem. For any path y : 7 — > R x R"^ 
withy(O) = (i,i/o)/ define /2y(s) := w{y{s)) and (py{s) := (p{y{s)). Since w is nonnegative, continuous, 
smooth on Of, and extends smoothly from Of to a function on R"^ for all t < 6, h'^^ :- lims^o+ ^y(s) 
and h'y_ := limg^o- ^y(s) exist. Simlarly, h'^^ and hy_ exist. Therefore, since w - cp is maximized 
at yo, h'y^ < (/)j,(0) < h'y_. Now, since w is not twice differentiable at yo, either there exists y such 
that hv is not differentiable at or if is differentiable at yo, but not twice differentiable. 

Case 1: y{s) is a path through yo for which hy is not differentiable. 
Then h'_ < h'^ and there exists no such (p. 

Case 2: w is differentiable at xq but not twice differentiable. 

Then, for all y through yo, (Py{0) = h'^XO) and (p"(0) > max(//^';+,?z^';_). Now, let y, be the 
coordinate paths through xq with w{yi{t)) > for < f < 6. Then, since on if > 0, if is a 
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subsolution of the linearized problem ( |4.1| |, we have 

(hdicj) - P{x, hD)cl)){xo) - hdfW - h^A<p - ih{Vp, D)w - jiw + V{x)w 



< hdfW - 



Llim h''. 

V i 



ih{yp, D)w - iiw + V{x)w < 0. 



Thus, we have that w is a subsolution on dOt- Putting this together with the arguments above, 
we have that w is a viscosity subsolution for ( |4.1| on i < 6. 

Now, by an ada ptati on of the maximum principle found in jH Section 3] to parabolic equations, 
any solution, v to (4.1 1 with initial data vq > wq will have v > w for t < 6. But, since v > 0, > 
and hence the solution u to with initial data vq will have u > v >w for t < 6. Now, since for 
f > 2, w{x, t) >1 on (pt{B{xo,a)), we have the result. 

We now construct a function vq satsifying the assumptions, (4.2|,( |4.3| l, and ( |4.4| |. Let Vi be the 
ground state solution of the Dirichlet Laplacian on B(0, 1) c R*^ i.e. 

-Avi = Avi on B(0, 1) Vild^^^^^^ = 0. 

Then, vi extends smoothly off of B(0, 1) and has i^i > inside B(0, 1) (see for instance IZl Section 
6.5]). 

Let X e C°°(B(0, 1)), < a: < 1 with a: = 1 on B(0, 1) \ B(0, 1 - e) and supp x c B(0, 1) \ B(0, 1 - 2e). 
Then, define V2 :— Mvi + [x{x)] {\x\^ - 1) where M is large enough so that V2 > 0. There exists such 
an M since z;i > in B(0, 1) and A > imply that -Avi = Avi > and hence, by Hopf's Lemma, 
dyVi < on (9B(0, 1) where v is the outward normal vector to dB{0, 1) (see for instance \Zt Section 
6.4]). 

For |x| < 1 - e, there exists C > such that, 

- Av2 = AMvi - {Axix)] {\x\^ - 1) - A{dx,x) - 2xd < AMvi + C. 
For |x| < 1 - e, Vi > 6. Thus, by increasing M if necessary, we obtain |S > such that 

-Av2 < AMvi + C < 2Av2 - /3. 

Now, for 1 - e < < 1, 

-Av2 = AMvi - 2d = A{V2 - \x\^ + l)-2d< Av2 + A(2e - e^) - 2d. 
Thus, for e > small enough, there exists fi > such that 

-AV2 < AV2 - ji- 

Finally, 3 a eR,xo e M!^, and Ci, C2 > constants so that 



Vo = 

satisfies the conditions on wq. 



^2i'V2{a ^{x - xq)) X e B{xo,a) 
else 



□ 
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Remark 1. If a shorter time is desired, one may use initial data of 0{h") to obtain the same result 
in time 0{h\ \ogh\). 

Remark 2. Notice that to obtain a growing subsolution it was critical that [i > 0. This corresponds 
precisely with the movement of the pseudospectrum of P{x, hD) into the right half plane. 

Now, we will demonstrate finite time blow-up using the fact that in time 0(1) the solution to 
is > 1 on an open region. The proof of theorem [l] follows 




Figure 5. We show simulations of solutions to the equation hut - Pi (x, hD)ii with 
h - 1/193. The solution using a quasimode mq with eigenvalue z - ^ and error 
0{h^) as initial data is shown in the blue solid linse. The red dotted lines show 
u{x, t) = uq{x) exp (zt/h). We see that the solution to the linearized problem ( |4.1| | 
with quasimode initial data closesly approximates the exponential until t « 0.3 

Proof. Let uo{x) and ti be the initial data and time found in Lemma |3] with {a,xo,6) such that 
<pt is defined on B{xo,a), (pt{B{xQ,a)) c V"^[0, |] for t e [0,6], and h < 6. Then, u{x,ti) > 1 on 
(pt{B {xo,a)). 

Now, let O € C^(]R) be a smooth bump function with 0(i/) = 1 on < 1, < O < 1, supp 
O c (-2, 2), and O" < CO^^^ (one such function is given by e'^^^ for e > x > 0). Define x '■^'^ 
hyxiy) :=o(2fl-i|y|). 



NONLINEAR INSTABILITY IN A SEMICLASSICAL PROBLEM 

Next, let y' = (pt{xo + y) and let 
Then, we have that 
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hvt - h^Av + fiv + v^ - 2h^{SIx, Vw) - h^Axu + (x - X^)"^ - ^(j/')^- 
Finally, define the operations, [/] and [/, g\ by 

[/] := f f{y)dy and {f,g] - f {f{ylg{y))dy. 

J B(0,fl) J B(0,fl) 

Then, we have that 

h{v]t = h^[Av] + ^[v] + [v^] - [Ax, u] - 2h^ [Vx, Vw] + [x- x\ - [V(y'). 
(4.5) > ii{v] + [v'] + [Ax, u] + [x- x\ - [Viy'l v] 

Here, (4.5 1 follows from integration by parts, the fact that = at |i/| = « and that 

f Av=^ \ Vv-v = Q. 

J B(0,fl) « J,9B(0,fl) 

We will later need that {v^} > [v]^. To see this use Holder's inequality with / = wc^{a)~'i, 
g - ,p - 3, and q - jto obtain 

- r„i3 



JB(0,a) \jB{0,a) / 



Next, we will estimate [Ax, u\. Observe that 



{Ax,u\ = 



^ — 0'(2fl-V) +4fl-V 



(4.6) 



(4.7) 



< C r (0^/2 + oi/3\ r u{r,(p)dSdr 

Jo J5B(0,fl) 



"(2fl-v)) r 

/ Jdl 



dB{0,r) 



u{r, (p)dSdr 



where C and C do not depend on h. Here ( |4.6| follows from the fact that for any function O > 0, 
O' < 0^/2 and that O' = near r = 0. follows from < O < 1. 

Now, we have 

h[v]t = fi[v] + + [Ax, u] + [x- x\ - [V{y'), v] 

> + [v^] - Oih^) + [(1 - Oih^))x - x\ - [Viy'), v] 

(4.8) > n[v] + [v^]-0{h^)-0{h^)[v^]-[V{y'),v] 

(4.9) > ^[v] + il-Oih^))[vf-Oih^)-[Viy'),v] 



Here, (4.8 1 follows from the fact that x ^rid (|4.9| follows for < 1 since [c ] > [c] . 
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Now, ont <6, we have V{y') < |. Thus, for < f < 6, 

h[v]t>^[v] + il-0{h^))[vf-0{h^). 

We have that [c](fi) > 1/4 and /j > 0. Therefore there exists y > such that, for h small enough 

and ti < t < ti + y, 

h[v]t>^[v] + ^[vf. 

But, the solution to this equation with initial data [v]{0) > 1/4 blows up in time t2 = 0{h). Hence, 
so long as ti + t2 < min(6, ti + y) and h is small enough, [v] blows up in time ti + ti- Observe that 



since t-[ < 6, < ti + t2 = ti + 0{h) < min(6, ti + y) for h small enough. Thus, the solution to 1.1 



blows up in time 0(1). □ 
Remark. A similar result holds for polynomially small data with blow up in time 0{h\ log h\). 



5. Numerical Simulations 

We expect that the instability of ( |1.1| is related to the presence of pseudospectrum in the right 
half plane. In fact, using numerical simulations for ( |4.1| based on code from |9 1 with P as in ( |1.5| , 
(see Figure |4]| we are able to demonstrate that the the solution with a quasimode for a positive 
eigenvalue as initial data closely approximates an exponential. Based on these results we expect 
that a proof of blow-up using quasimodes will allow the results of Theorem [l] to be extended to 
complex energies and wider classes of operators. 

All simulations were run in 1 dimension with = | • Unless otherwise stated, all quasimodes 
are constructed with errors of 0{h^). 
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